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In nonuniform Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) superconductors, both the gauge sym- 
metry and the continuous translational symmetry of the normal state are spontaneously broken. 
This leads to additional bosonic excitations, or Goldstone modes, corresponding to the deforma- 
tions of the order parameter amplitude modulation in real space. We derive general expressions for 
the energy of the phase and elastic Goldstone modes. As an example, the superfluid density and 
the elastic modulus of a one-dimensional LOFF superconductor are calculated at low temperatures. 
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I. INTRODUCTION 

In most superconductors, the Cooper pairing is suppressed by the coupling of the orbital motion of electrons 
with a magnetic field. 1 In some cases, however, the orbital suppression is ineffective and the dominant mechanism 
of the magnetic pair breaking is spin-related, caused by either the Zeeman coupling with the external field, or 
the exchange interaction with the localized spins in a magnetic crystal. The superconductivity is then said to be 
paramagnetically limited. It was shown by Larkin and Ovchinnikov 2 (LO) and Fulde and Ferrell^ (FF) that the 
competition between the paramagnetic pair breaking and the condensation energy results in the formation of a 
peculiar nonuniform superconducting state - the LOFF state - in which the Cooper pairs have nonzero center-of- 
mass momenta. In the simplest cases, the order parameter is described by either a single plane wave: A(r) = Aoe 1 ''*' 
(FF state), or a superposition of two plane waves: A(r) = AoCosQr (LO state). In general, minimization of the free 
energy yields more complicated structures. 

The experimental detection of the LOFF state requires both the orbital effects and disorder to be sufficiently 
weak. This can be achieved in the superconductors with a short coherence length and heavy effective mass of 
quasiparticles, which increases the orbital upper critical field H C 2, and can be particularly favored in the systems 
with low-dimensional geometry, either in real space (superconductivity in films or at interfaces, in a parallel field) or 
in momentum space (quasi-one-dimcnsional and quasi-two-dimensional band structures). For a review of the LOFF 
state in superconductors, see Ref. 0. A different route to the experimental realization of the LOFF state is offered 
by cold atomic Fermi gases. In these systems the role of the Zeeman field is played by the mismatch between the 
Fermi surfaces, which is controlled by the difference between the numbers of atoms in two different hyperfine states, 
see Ref. HI for a review. Further afield, states similar to the LOFF state have been discussed recently in high-energy 
physics, in the context of "color superconducting" quark matter £ While most theoretical studies have focused on the 
equilibrium properties of the LOFF systems, in particular, the phase diagram and the most stable order parameter, 
a number of recent works have looked at deviations from a perfectly periodic order, including thermal and quantum 
fluctuations^— and topological defects, such as vortices and dislocations^ : 11 ' 12 

In addition to the U(l) phase rotation symmetry, the LOFF state can break the continuous translational symmetry, 
which leads to a richer spectrum of low-energy bosonic excitations, or Goldstone modes. Since neither global phase 
rotations nor uniform translations of the order parameter cost any energy, the Goldstone modes in the LOFF state 
are associated with slow gradients of the order parameter phase, as well as with weakly nonuniform displacements of 
the order parameter amplitude modulation. In the lowest order, the energy is quadratic in both the phase and the 
displacement gradients. For instance, one can expect that the energy of long-wavelength static phase fluctuations has 
the usual form p s w 2 /2, where the superfluid velocity v s is proportional to the gradient of the fluctuating phase, and 
p s is the superfluid density, which characterizes the phase stiffness of the superconducting state. 

In this paper, we present a microscopic derivation of the free energy of the Goldstone modes in a LOFF super- 
conductor. Recently, a similar problem has been considered using the Ginzburg-Landau expansion near the critical 
temperature in the isotropic case, which is applicable to either a Fermi gas with a population imbalance^ or the 
quark matter^ In Ref. [l5|, the elastic moduli of a LOFF vortex lattice were calculated. In contrast to these works, 
which used the free energy expansion in powers of the order parameter, we employ a transformation of the slow de- 
formations of the LOFF order parameter into small perturbations in the effective bosonic action. General expressions 
for the effective action of the static Goldstone modes, both in the FF and LO states, are derived in Sec. |TTJ In Sec. 
IIII1 we illustrate how the formalism works using as an example the LO phase in a one-dimensional superconductor at 
low temperatures, for which the order parameter and the spectrum of excitations in the mean-field state are known 
exactly. Throughout the paper we use the units in which h = ks = 1. 



2 



II. EFFECTIVE BOSONIC ACTION IN THE LOFF STATE 



We consider a spin-singlet superconductor without disorder, in an external magnetic field B. The orbital pair 
breaking is neglected, so that the superconductivity is affected by the field only through the Zeeman splitting of the 
electron bands. The Hamiltonian is given by 



H 



dripi(r)(£5 a p + ha 3ta p)ipp(r) - V f dri/j\(r)ipl(r)ip±(r)ip^(r). (1) 



In the first term, £ = £(fc) is the effective band Hamiltonian (which includes the chemical potential /x), k — — iV, 
a,/3 =t,4- is the spin projection on the quantization axis along B, h — (g/2)fj,sB is the Zeeman field (the electron 
charge is equal to — e), /jb is the Bohr magneton, g is the Lande factor, and 03 is the Pauli matrix. The electron 
wavefunctions are assumed to satisfy the periodic boundary conditions in a cubic box of side L and volume V = L 3 . 
The second term in Eq. (JTJ) is the Bardeen-Cooper-Schrieffer (BCS) pairing interaction with the coupling constant 
V > 0. The Hamiltonian (JTJ can also be applied to a ferromagnetic superconductor in zero applied field, in which 
the electron bands are split due to the exchange interaction with the magnetically ordered localized spins, or to an 
imbalanced atomic Fermi gas. 

We derive the free energy of our superconductor starting with the standard expression for the partition function in 
terms of the Grassmann functional integral: Z = J T>^) a T>^) a e~ , see, e.g. Ref. [l6f. Here ip a (r,r) and ^j a (r,T) are 

fermionic fields, S = Jq dr [J dr -ipadripa + H{t)\ is the action associated with the Hamiltonian ([1]), and f} = 1/T. 
Decoupling the interaction term by the Hubbard-Stratonovich transformation, the partition function can be written as 
a functional integral over a complex bosonic field A(r, r) which has the meaning of a fluctuating superconducting order 
parameter: Z = J VAVA* e~ s " ff ^ A \ where S e ff is the effective bosonic action. Focusing on the static fluctuations, 
the effective action takes the following form: S e ff — )3J-[A], where 

J-=-TVTrlnCr 1 + -^ [d 3 r\A(r)\ 2 (2) 



n 

is the free energy, Q~ x = icu n — 7-{, u) n = (2n + 1)ttT is the fermionic Matsubara frequency, and 

U=(it k f*} (3) 



A*(r) -i + h 

is the Bogoliubov-de Gennes (BdG) Hamiltonian. The trace in the first term in Eq. (|2|) is taken in both the 
coordinate and the electron-hole spaces and can be formally written as follows: Tr lnty 1 = ^ a ln(za; n — E a ), where 
E a are the eigenvalues of %, labelled by quantum numbers a. The sum over a formally diverges and should be 
regularized by subtracting the corresponding expression in the normal state (in the spirit of the BCS model, we 
assume that only a finite number of the electron states in a narrow energy interval near the Fermi energy are affected 
by superconductivity). Since we are interested in the energy of fluctuations, which vanishes in the normal state, we 
will not write the regularizing terms explicitly. 

The mean-field order parameter, denoted by Ao(r), corresponds to a saddle point of the effective action, which 
is found from the equation ST/SA* = 0. The solution is sensitive to the electronic band structure and the pairing 
symmetry. The problem of finding the most stable LOFF state at all T and h has not been solved even in the simplest 
cases, apart from an exact solution in one dimension, see Refs. [T3 and [l8l . In a fully isotropic three-dimensional 
system, it was argued in Ref. [l9| that the phase transition between the normal and the LOFF states is always first 
order, and that the order parameter is represented by a sum of one, two, or three cosines, as the temperature is 
lowered. In this work, we consider only the phases with a one-dimensional (ID) periodicity, see Appendix [XJ the FF 
phase with 

A ,FF(r) = A e lQz , (4) 

and the LO phase with 

A .Lo(r) = A /(z), (5) 

where f(z) is a real periodic function with the period d. In the vicinity of a second-order normal metal-LOFF 
superconductor phase transition, one can put f(z) — cosQz, where Q = 2n/d. Away from the critical temperature, 
nonlinear effects add higher Fourier harmonics to the order parameter: 

00 

/(*)= E fp eVQz > /-* = /;• (6) 

p— — 00 
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The sum here does not include the p = harmonic, so that the spatial average of the order parameter is zero. Thus, 
the nonuniform superconducting state resembles a periodic array of domain walls separating the regions where the 
order parameter is almost uniform.— 



A. FF phase 



The mean-field solutions described above are not unique. While one can uniformly rotate the overall phase of the 
order parameter without any energy penalty, a nonuniform phase rotation corresponds to a Goldstone mode, whose 
energy is small in the long- wavelength limit. In contrast, even uniform deviations of the order parameter amplitude Ao 
from its mean-field value cost energy and can be neglected. Thus, the low-energy fluctuations of the order parameter 
in the FF phase are described by the following expression: 



A FF (r) = A Q e w( ^e lQz , 



(7) 



where varies slowly on the scale of the LOFF period d. The fluctuating phase can be removed from the off-diagonal 

elements of the BdG Hamiltonian by the unitary transformation T~L — » fi = cxp(—i8a 3 /2)'Hexp(i9d' 3 /2), and we 
obtain: 



H = 



£(tt+) + h A e^ z 
Aoe" jQz -£(*_) + ft 



(8) 



where -k± — k ± mv s , v s = V6>/2m is the superfluid velocity, and m is the electron mass. 

For slow fluctuations, one can use a constant superfluid velocity and expand the BdG Hamiltonian 

of v s . For the lowest two orders in the perturbation series we obtain: % = T-Lq ff + 5%, where 



FF 



i + h A„e iQz 
A e-^ z -£ + h 



is the mean-field Hamiltonian, 



m powers 



(9) 



SU = mv s .iVi(k)cr Q + -m 2 v s ^v s jm il 1 (k)a 3 , 



(10) 



v = d^/dk is the quasiparticle velocity, and = 3 2 £/ dkidkj is the inverse tensor of effective masses. 
Using Eq. the phase fluctuation contribution to the free energy can be written as follows: 



5F FF = -T^[Tr ln^o" 1 - 8U) - Tr ln^ 1 ] = 



tt (g sH) + -ti (g sng SH) 



where C?o = (ik> n — ?^o,ff) _1 is the 2x2 matrix Green's function in the undeformed FF state. Keeping only the 
quadratic terms in v s (the first order terms vanish at the saddle point), we obtain the energy of the phase Goldstone 
modes per unit volume: 



FF 



V 



1 

'^Ps.ij'^s.i'^s.j i 



where 



£ [Tr {vAvjQo) + Tr (m^ 3 6o)] 



(11) 



fl2) 



is the superfluid mass density tensor. In the normal state, the Green's function is electron-hole diagonal, 
translationally-invariant, and satisfies the identity 



dg (fc,q;„ 
dk 



v(k)a 3 g^(k,uj n ). 



(13) 



Using ni- 1 = dvi/dkj and integrating the second term in Eq. (I12[) by parts, one can verify that in the normal state 
the superfluid density vanishes. 
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B. LO phase 

In addition to the phase rotation symmetry, the LOFF order parameter can also break continuous translational 
symmetry, leading to the existence of additional Goldstone modes, see Appendix El In the LO phase, these additional 
modes correspond to weakly nonuniform deformations of the amplitude modulation. 

Suppose a Cooper pair that was at a point r' in the undeformed state is found at a point r after the deformation 
characterized by the displacement ix(r') = r — r' (we assume that the relation between r and r' is unique and 
invertible). As a result, the pair wavefunction Ao(r') is transported to the new location, leading to the order 
parameter at the point r being transformed from Ao(r) into A(r) = Ao(r'(r)). For the LO phase with a ID 
periodicity, described by Eq. ([5]), the deformed order parameter has the following form: 

A LO (r) = A e ie ^f(z\r)), (14) 

where the relation between r and r' is defined by the expressions x = x' , y = y' , z = z' + u{x' ,y' ', z'), and both 
6 and u vary slowly compared with d. Due to the periodicity conditions for the wavefunctions, the region 1Z in 
three-dimensional Euclidean space occupied by the superconductor is homeomorphic to a three-dimensional torus. 
The "Eulerian" coordinates r and the "Lagrangian" coordinates r' (Ref. |2l|) correspond to different parametrizations 
of the points in 1Z. Since the order parameter also has to be periodic, we assume that u[x' + i, y', z') = u(x', y', z'), 
etc, therefore, x'(x + L,y,z) — x'(x,y,z) — L, etc. Since neither a global phase rotation, described by a constant 
0, nor a uniform translation, described by a constant u, cost any energy, the free energy T can only depend on the 
gradients of 9 and u. 

It is convenient to transform the BdG Hamiltonian of the LO phase into the coordinate system in which the order 
parameter is not deformed. This can be achieved by changing coordinates from r to r' , which affects the metric of 
the Euclidean spaced In the Lagrangian coordinates, Eq. © takes the form 

y, = ( i' + h A e i9 ^f(z') 
\A e- i6 ^f(z') -i' + h 

while the components of the metric tensor become 

4, + Sk3-^-j [S kj + fe^T-J , (16) 

where we introduced the notations (a; 1 , x 2 , x 3 ) = (x,y,z) and (a;'' 1 , a; 7,2 , a; 7,3 ) = (x',y',z'), and the repeated indices 
are summed over. In order to find how the band Hamiltonian is transformed, we assume that it can be represented 
as a series expansion: 

£ = y y A 1 N ki 1 ...ki N — fx, 

even iV 

where are determined by the crystal symmetry and completely symmetric with respect to the permutation of 

«i, ijy = 1,2, 3. Under the change of coordinates r — > r', this is transformed into 

^' y.-^ ■'<;'<■■■'::><, /'• (17) 

N 

where k[ — — zVJ and e^(r') = dx'^ jdx % . For instance, if the band dispersion can be treated in the effective mass 
approximation, i.e. £(fc) = fc 2 /2m* — fi, then |' = — A LB /2m* — [i, where A LB = .9 _1/ ' 2 V-g 1 / 2 g lJ ; Vj- is the Laplace- 
Beltrami operator (Ref. 22), g = (1 + du/dz') 2 is the determinant of the metric tensor, and are the components 
of the inverse of the metric tensor. 

The partition function can be written as follows: Z oc j VAVA* e^^^ 1+jr ' 2 \ where 

e -/»* = JvVVVe^^' 1 ^, (18) 

^ = (^,i/i^) T and ^ = (ifo,ipi) are the Nambu Grassmann fields, T% = V^ 1 J d 3 r |ALo( r )| 2 j and A^o is given by 
Eq. ([T4"]) . One can show that does not depend on the deformation. Indeed, using Eqs. (J5J) and ©, we obtain: 

= i / d 3 r'V5|A ,Lo(r')| 2 = ^ E ^vJvJv,- 

P1P2 
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Here I PlP2 = f d 3 r' y/g e^+^Qz' ~ V5 Pu - P2 , if the displacement u varies on the scale much greater than the LOFF 
period. Therefore, Ti = (VA(j/V)5^ \f P \ 2 , i-e. all the effects of the order parameter deformation are contained in 
./•i • 

Since the transformation from r to r' changes the scalar product in the functional space, one should be careful 
calculating the integral in Eq. (|18[1 . Let us introduce a complete and orthonormal set of wavefunctions in 1Z, using the 
Eulerian coordinates r. One can use, for instance, the eigenfunctions of £, i.e. the normalized plane waves Xfc( r ) = 
V~ 1 / 2 e lfc ' r , where k — (2n/L)(mi, 1712,1713), and nrii take integer values due to the periodicity conditions. We construct 
the basis Nambu spinors labelled by a — (k,s) (s = 1,2 is the electron-hole index), as follows: <ph,i = (Xk,0) T and 
<Pk,2 = (0, Xk) T , an d represent the Nambu fields in the form 

n a n a 

where c n ^ a and c n ^ a are elements of the Grassmann algebra. From Eq. f|18p we obtain: 

e-PK = f JJ JJdc^^eEn..*^ 1 ^-)^..^ = FJdetg- 1 . (19) 

n a n 

The matrix elements of the inverse Green's operator are given by 

GabM = f d 3 rtr^( r )( lWn _H)^( r)= f dV^tr<^( r ( r '))^„_HV(KO), 



where "tr" denotes a 2 x 2 matrix trace and H' is given by Eq. (fT5|) . Introducing the new Nambu wavefunctions 

<p a (r') = g^e-WWipaW)), (20) 

we obtain: 

Q'bM = J dVtr^ (r ')( iWn -H)(p b {v'), (21) 

where 

H M+ + h ) (22) 



and 



& = e ^/y/ 4 £V 1/4 e^ (23) 

Af 

The last expression follows from Eq. (fl7|) . 

The functions (|2ll|) form a complete and orthonormal set in TZ. Since the factors g 1 ^ 4 have been absorbed into the 

definitions of (p a , the integral in Eq. (|2ip is formally the same as for fermions with the Hamiltonian % moving in a flat 
space (with the metric tensor given by a unit matrix). Now one can drop the primes and write <p a ( r ) — Y]h fbi^Uba, 
where J7(, a = J d 3 r tr <p b (r)<p a (r) form a unitary matrix. Therefore, using Eq. p^|) . we obtain: detQ^ 1 = det[W(iuj n — 

H)U] = det(io; n — %), and 

T x = -T^2 Tr ln(iuj n - H). (24) 

n 

The advantage of using the transformed Hamiltonian H is that, in contrast to ~H, it contains only the gradients of the 
fluctuating fields, which are small and can be treated perturbatively. 

We consider only the limit of a slow deformation, when the gradient of the displacement can be set to a constant: 
Vit = £. Different components of e play different roles: while e XtV correspond to a uniform tilting, e z describes 
a uniform compression or stretching of the order parameter modulation. The effective band Hamiltonian (|17[) is 
transformed into £ = £(fe — ek z /(l + e z )), while Eq. (|2"3"|) takes the form £± = £(7r± — en± iZ /(l + e z )), where 
tv± — k ± (recall that the primes denoting the Lagrangian coordinates have been dropped). Although the 
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assumption of a constant strain is not consistent with the global periodicity conditions, it is legitimate if the scale 
of variation of u is much greater than d, the period of the LOFF structure. This "local limit" fails when d diverges, 
which happens, e.g., near the tricritical point on the T — h phase diagram, where the Ginzburg-Landau gradient term 
changes sign. Similarly to the FF case, we also assume that the superfluid velocity v s is uniform and small. 

For slow fluctuations, we expand the Hamiltonian (j2"2")l in powers of v s and e, and obtain: H = H lq + 8H, where 



is the mean- field part, while 



U °> L °- [ A„/(z) (25) 



5% = (lo 1 - e r LLi z )vi(k)a + ^tiiLbjm^ 1 (k)a- 3 , (26) 



with u) = mv s ao — ek z a 3l can be treated as a small perturbation. We keep only the lowest two orders in v s and e. 
The higher orders can be easily obtained from Eq. (f2U)l , if needed. 
The energy of the Goldstone modes follows from Eq. (|24|): 



where C?o — (*w n — Hq.lo) 1 the 2x2 matrix Green's function in the undeformed LO state. Retaining only the 
quadratic terms, we arrive at the following expression: 

Here the superfluid density tensor is given by 



p s ,ij "' J y 



gT^[Tr(v^{i J -0 o ) + Tr(mr. 1 ^o)] ) (28) 



which has the same form as in the FF state, see Eq. (TT21 . but with a different Green's function, 

Kij = — T^2[Tr (vik z a3QoVjk z o-3Go) + Tr (rh'j 1 kla 3 Q Q ) + 2S lz Tr (vjk z a 3 go)] (29) 

n 

is the tensor of elastic moduli, and 

Kij = ~m^T^2[TT (v i g Q v J k z a 3 g ) + Tr [m^k z Qo) + (5 iz Tr (vjQo)] ■ (30) 

n 

It is straightforward to show, using the identity (p~3|) . that all the stiffness coefficients, Eqs. (|28|) . (|29|) . and (|30l) . 
vanish in the normal state. 

If the order parameter has a center of inversion, i.e. A(— r) = A(r), then the Green's function satisfies 
^o(— ri, — r 2 ;o; n ) = ^o( r i, r 2', w n ), and Eq. (I5TH) yields — 0. Therefore, the Goldstone modes corresponding 
to the phase fluctuations and elastic deformations are decoupled. This can also be understood using the following 
symmetry argument. Assuming a general three-dimensional displacement u = (u 1 , it 2 , it 3 ), one can introduce the 
displacement gradient tensor — du 1 jdx^ (the usual strain tensor of the elasticity theory is the symmetric part of 
ij). For a uniform deformation of the LO phase, we have u 1 = u 2 = 0, u 3 = e\x l , where e = (ef , e^j ^3) does not 
transform like a vector. In particular, e remains invariant under inversion, while v s changes sign, therefore the free 
energy cannot contain quadratic terms mixing e and v s . 



III. LO PHASE IN ONE DIMENSION 



As an application of the general formalism developed in the previous section, we consider the LO phase in a 
superconductor with a ID band dispersion £(fc z ) = {k 2 z — k 2 F )/2m*, where ftp is the Fermi wavevector and m* is the 
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FIG. 1: The order parameter (solid line) and the quasiparticle wavefunctions (dashed line) in the ID LO phase. 



effective mass. In this case, the mean-field gap equations can be solved exactly, see Refs. [l7j and[l8|, using the formal 
similarity with the ID Peierls problem considered in Ref. [23l 

Let us briefly summarize the relevant properties of the exact solution. Although the critical temperature decreases 
with the Zeeman field h, the superconductivity is not completely suppressed even at strong fields. The order parameter 
is described by the following expression: 

A(z) = A 1 cd(-^-z,k 1 ), (31) 

where "cd" is a Jacobian elliptic function (Ref. HJ), vp = kp/m* is the Fermi velocity, and Ai and k\ are parameters 
that depend on T and h. The period of the LO structure is given by 

where 7 = 2\fk\j{\ + fei), and K(7) is the complete elliptic integral of the first kind. At low temperatures, the LO 
structure resembles a soliton lattice, sketched in Fig. [T] Note that our expression for the LO order parameter, Eq. 
pip , is shifted by a quarter-period compared to those in Refs. [I?] and HH to make explicit the inversion symmetry 
A(— z) — A(z). At zero temperature, one has Ai = \fk[Ao, where Ao is the BCS gap at T = h = 0, while the 
parameter 7 as a function of h is found from the equation E(7)/7 = 7rft,/2Ao, where E(7) is the complete elliptic 
integral of the second kind. 

The energy density of the Goldstone modes, Eq. (|27p. takes the following form: 

with the coefficients in the first and second terms characterizing the stiffness of the LO phase against the phase 
fluctuations and the deformations of the amplitude modulation, respectively. We have p s = p\ + pi , where 

ft = 4 T E7 Tr ^o), P2 = ^Tj2T r b(kMJo), (34) 

n n 

and K = Ki + K2, where 

Ki = -^rrk^o), k 2 = J^TY^jTr(k%g k%g ) (35) 



(L is the system length). 

Since the order parameter varies slowly compared to the Fermi wavelength hp 1 , one can use the quasiclassical, or 
Andreev, approximation. The wavefunctions can be represented in the form e lskFZ ip(z), where s = ± is the direction 
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FIG. 2: The quasiparticle excitation spectrum in the ID LO phase, shown in the extended zone scheme (from Ref. [23). The 
Zeeman field h is located inside the energy gap. 



index, which labels the roots of the equation £(k z ) 
the Andreev Hamiltonian: 



0, and the slowly varying factors ip(z) are the eigenfunctions of 



H s = —isvFO'3- — h A(z)a\. 
dz 



(36) 



Due to the periodicity of A(z), the wavefunctions satisfy ip(z + d) = e lq tp(z) and are characterized by the wavevector 
q. At given s and q, the spectrum consists of two branches, labelled by the branch index v — 1,2, with E s q 2 ~ ~E qll 
see Fig. [2] There are two energy gaps, located at q = ±w/d, with the gap edges given by ±e + and ±e_, where 
e_ = Ao^/l/7 2 — 1 and e + = A0/7 (at T = 0). It is easy to show that the Zeeman field is located inside the gap, i.e. 
e_ < h < e + , and that e + — h < h — e_. The matrix Green's function can be written as follows: 



Go(z, 







isk F (z-z') ' 



-Ei-h' 



(37) 



The quasiclassical approximation cannot be directly applied to pi, see Eq. (|34[) . because of the contributions of 
the states far from the Fermi energy. To get around this difficulty, we subtract from and add to pi the corresponding 
expression in the normal state: pi = L' 1 (m 2 /m*)T ~^2 n Tr (a^Go.N) + Spi, where Q~^(k z , uj n ) — iui n — £,{k z )a 3 — h, 
while 



Spi = — -TV \ [Tr (& 3 Go) - Tr (<7 3 ^o,7v) 



can be calculated quasiclassically. If the asymmetry of the quasiparticle density of states near the Fermi energy is 
neglected, then Spi = 0. The remaining expression is entirely determined by the normal-state properties and yields 
pi = (m 2 /m*)n, where n is the concentration of electrons. 

The second contribution to the superfluid density, pi, comes from the states in the vicinity of the Fermi energy and 
can therefore be calculated in the quasiclassical approximation. Substituting the Green's function (|3T|) in Eq. (|34p . 
differentiating only the fast exponentials, neglecting the integrals of functions that oscillate on the scale of kp, and 
using the orthonormality of the Andreev eigenstates, we obtain: 



P2 



2 2 
-TO Vp 



dEN(E) 



df(E + h) 



dE 



(38) 
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where f(E) = (e f3E + is the Fermi distribution function, and N(E) = L^ 1 J2 S qu — Eq, v ) tne density of 
states of the Bogoliubov quasiparticles. At zero temperature, we have P2 — —m 2 v F N(—h) — 0, because —h is located 
inside the gap, see Fig. [2] At low temperatures, the corrections are exponentially small: piiT) oc e~ Eg / T , where 
E g {h) = e + — h = Ao[l — 2E(7)/tt]/7. Putting all the pieces together, we find 

Psi T) = ^n-0(e- E */ T ). (39) 

Thus the superfluid density at T = is the same as in a uniform BCS superconductor. At nonzero temperatures, it 
is reduced due to the thermally excited Bogoliubov quasiparticles, characterized by the energy gap E g . 

Let us now turn to the calculation of the elastic modulus, see Eq. (|3"5j) . As in the case of pi, the quasiclassical 
approximation cannot be applied directly to K\, but one can write K± = (3/m* L)T^2 n Tr (fc 2 <73<5o.jv) + SKi, where 
5K, = (3k 2 F /m 2 )S Pl = 0. Therefore, 

3 1 9k 3 

m *-~ / L * — ' nm 

n k z 

at zero temperature. Here we neglected the corrections of the order of h/ep, where cf — k F /2m*. The quasiclassical 
expression for K2 reads 

j4_ 1 l2 m qi ^ + h)-f(E^ 2 + h) 



J qi,Ui -"92, ^2 



where the matrix element is given by Ig lVl q2V2 — Jq dz ' i Pq' 1 \v 1 { z )G'3' l l>q 2 V J\ Z )- The interband contribution is dominated 
by the "vertical" transitions near the gap edges at q = ±7r/d, where the matrix element vanishes for symmetry reasons. 
Keeping only the intraband terms in Eq. (|4ip . we obtain: 



777, ' J 



'■' ' - Of{E + hy (42) 



dE 



where N{E) = L 1 Yl S qu \Iqv,qu?${E — Eq,v)- At zero temperature, we have K 2 oc N(—h) — 0, while the low- 
temperature corrections are exponentially small and proportional to e~ Ea ^ T . Thus, the elastic modulus of the ID LO 
state at low temperatures has the following form: 

Ob 3 

K (T) = fh^^O{e- E ^ T ). (43) 
7rm* 

Similarly to the superfluid density, it is reduced by the thermally activated Bogoliubov quasiparticles, until both p s 
and K vanish at the critical temperature T c {h). 

To conclude this section we note that the nonvanishing superfluid density p s implies that a uniform supercurrent 
can flow across the zeros of the order parameter amplitude. This can be understood using a simple phenomcnological 
argument based on the Ginzburg-Landau theory for the LOFF state. Near the critical temperature, the lowest order 
terms in the gradient energy can be written as follows: J- gra d = J dz (K2\D Z A\ 2 + A4|D 2 A| 2 ), where D z = — iV z , and 
K2 < 0, K4 > 0, in order for the superconducting instability with a finite wavevector to occur^ The supercurrent 
can be calculated in the usual way: j sz = —c(6F gra d/5A z ), after making the replacement D z — > —iS/ z + (2e/c)A z in 
■Fgrad, where A z is the vector potential and e is the absolute value of the electron charge. We obtain: 

j s , z = -4eK 2 Im(A*V z A) - 4eA 4 Im [(V Z A*)V 2 A - A*V^A] . (44) 

Using the amplitude-phase representation of the order parameter: A(z) = |A(z)|e* e( - z ), and keeping only the lowest 
order terms in the superfluid velocity, we obtain: 

A 2 |A| 2 + 2A 4 ((V Z |A|) 2 -2|A|V 2 |A| 



Thus, the supercurrent is not just proportional to |A| 2 and, therefore, does not vanish at the points where the order 
parameter amplitude has zeros. 
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IV. CONCLUSIONS 

We have derived general expressions for the energy of the Goldstone modes in the LOFF phases with one-dimensional 
periodicity, see Eq. (TTTj) for the FF phase and Eq. (|2"T)) for the LO phase. While there is only one Goldstone mode in 
the FF state, corresponding to the fluctuations of the order parameter phase, there are additional, "elastic" , modes 
in the LO state, corresponding to the compression and tilting of the order parameter amplitude modulation. Our 
approach, which is based on the transformation of slow deformations of the order parameter into small corrections 
to the Bogoliubov-de Gennes Hamiltonian, allows one to treat the fluctuation effects perturbatively. It can be easily 
generalized to more complicated LOFF structures. 

Using the exact solution for the excitation spectrum in a one-dimensional LO phase, we calculated the stiffness 
coefficients for the phase and elastic Goldstone modes. At T = the phase stiffness (the superfiuid density) has 
the same form as in a uniform superconductor, while at low temperatures the corrections are exponential, with a 
field-dependent energy gap. 

This work was supported by a Discovery Grant from the Natural Sciences and Engineering Research Council of 
Canada. 

Appendix A: Counting Goldstone modes in the LOFF state 

In the vicinity of the critical temperature T c (h), the mean-field free energy of a LOFF superconductor can be 
expanded in powers of the order parameter as follows: 

T = A{q)A* q A q + B(q ll q 2 ,q 3 ,q i )A* qi A* q2 (Al) 

q qi,2,3,4 

Since the full momentum dependence of A and B is retained, Eq. (|A1|) can be used close to the second-order phase 
transition line at all h. The coefficient A(q) changes sign at the critical temperature of the superconducting instability 
with the wave vector q. The explicit form of the functions A and B depends on the microscopic details. 

The order parameter in the LOFF state close to the critical temperature can be represented in the following form: 

Nq 

A(r) = ^A i e < ^, (A2) 

i=l 

where Aj = |A$|e"" are complex coefficients, Qi are the positions of the degenerate minima of A(q), and Nq is the 
number of the minima. In order to find the order parameter components A,, which determine the spatial structure of 
the LOFF phase, one needs to evaluate the free energy (|A1[) for the order parameter (|A2I) . We obtain: F /V = F 2 +F 4l 
where F 2 = A n J2i |Ai| 2 , with A — A{Q{) (same for all i), and 

F± = J2 B tjkl A*A*A k A ly (A3) 

ijkl 

with Bijki — B(Qi, Qj, Q k , Qi)5Q i +Q jt Q k +Q l ■ In addition to the momentum conservation and the conditions Hyfez = 
Bjiki = Bijik, the coefficients B^i must also satisfy certain symmetry-imposed constraints, which can be obtained 
from the requirement that the order parameter transforms like a scalar function under an arbitrary operation g from 
the point group G of the crystal, i.e. A(r) — > A(g^ 1 r). Inserting here the expansion (|A2[) and taking into account 
that the set of Qj's is invariant under all operations from G, we see that the action of g on the set of Aj's is equivalent 
to a permutation P(g). Therefore, Eq. (IA3|) must remain invariant under P(g) for all g. 

Let us illustrate the above statements using as an example a tetragonal crystal with G = D^. The point group 
has 16 elements, therefore there can be as many as Nq — 16 degenerate minima of A(q). The general case is clearly 
untreatable, so we consider just two simplest cases: Nq — 2, in which the deepest minima of A are located at 
Qi t 2 = ±Q, where Q = Qz; and Nq = 4, in which the minima occupy four high-symmetry points in the basal plane: 
Q 12 = ±Qx, Q 3 ,4 = ±Qy- 

Nq = 2. The order parameter has the form A(r) = Aie l ® z + A 2 e~ lC ^ z . The free energy must be invariant under the 
action of the generators of D^, i.e. rotations C^ z and C 2x , and also inversion /. Due to the momentum conservation, 
the nonzero coefficients in Eq. (|A3j) are Bun = B(Q, Q, Q, Q), B 12V2 = B 2112 = B 1221 = B2121 = B(Q, —Q, Q, —Q), 
and B2222 = B(—Q, —Q, —Q, — Q). Since the inversion operation interchanges Ai and A2, we have Bun = B 2222 , 
and Eq. (|A3[) takes the following form: 

F 4 = /MIA1I 4 + |A 2 | 4 ) +/? 2 |A 1 | 2 |A 2 | 2 , (A4) 
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where fa = Bun and fa = 4i?i2i2- The free energy does not depend on the phases of A1.2, therefore there might 
be up to two Goldstone modes in the LOFF state, corresponding to the fluctuations of either the overall phase of 
the superconducting order parameter or the relative phase of Ai,2- Expression (|A4[) is positive definite if fa > 0, 
fa > — 2/8].. If fa < 2/3], then the minimum energy is achieved for |Ai| = IA2I = Ao, which corresponds to the 
LO phase with A(r) = Aoe i6 cos(Qz + ip). The relative phase of A]^ describes a uniform translation of the order 
parameter. On the other hand, if fa > 2/3], then IA2I = or |Ai| = 0, corresponding to one of the two degenerate 
FF phases, A(r) = A e ie e i( ^ z , or A(r) = A e' e e- lC ^ z . 

N Q = 4. The order parameter has the form A(r) = A ie lQx + A 2 e~ lQx + A 3 e lQy + A A e~ iQv . Under C iz , we have 
A] 2,3,4 — > A3 4 2.1, while under C 2x and I, A] «-» A 2 and A 3 f-> A4. The nonzero coefficients in Eq. (|A3|) are as 

follows: i?nn = -B2222 = B 3333 — -B4444 = fa, i?1212 = #3434 = fa/A, -B1313 = i?1414 = ^2323 = #2424 = fal^i and 

£1234 = ^3412 = At/ 4 , therefore, 

F 4 = /3 1 (|A 1 | 4 + |A2| 4 + |A 3 | 4 + |A4| 4 ) + /32(|A 1 | 2 |A2| 2 + |A 3 | 2 |A4| 2 ) 

+fa(\A 1 \ 2 + |A 2 | 2 )(|A 3 | 2 + |A 4 | 2 ) + &(A; A*A 3 A 4 + A]A 2 A*A:). (A5) 

The last term leads to a phase locking of the order parameter components: at the energy minimum the phases of Aj 
are subject to the constraint (f>\ + (f> 2 — 4> 3 — 4>4 — or n, depending on the sign of fa. Therefore, there can be up to 
three Goldstone modes, one of which corresponds to the fluctuations of the overall phase of A(r). 
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